Selection rules for Single-Chain-Magnet behavior in non-collinear Ising systems 
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The magnetic behavior of molecular Single-Chain Magnets is investigated in the framework of a 
one-dimensional Ising model with single spin-flip Glauber dynamics. Opportune modiflcations to 
the original theory are required in order to account for reciprocal non-collinearity of local anisotropy 
axes and the crystallographic (laboratory) frame. The extension of Glauber's theory to the case of 
a coUinear Ising ferrimagnetic chain is also discussed. Within this formalism, both the dynamics 
of magnetization reversal in zero field and the response of the system to a weak magnetic field, 
oscillating in time, are studied. Depending on the geometry, selection rules are found for the 
occurrence of slow relaxation of the magnetization at low temperatures, as well as for resonant 
behavior of the a.c. susceptibility as a function of temperature at low frequencies. The present 
theory applies successfully to some real systems, namely Mn-, Dy-, and Co-based molecular magnetic 
chains, showing that Single-Chain-Magnet behavior is not only a feature of coUinear ferro- and 
ferrimagnetic, but also of canted antiferromagnetic chains. 

PACS numbers: 75.10.-b, 75.10.Pq, 75.50.Xx, 75.60.Jk 



I. INTRODUCTION 



Slow dynamics of the magnetization reversal is a cru- 
cial requirement for potential applications of Single- 
Chain Magnets (SCM's)i>^!^, and nanowires in general, 
in magnetic- memory manufacture. For nanowires with 
a biaxial anisotropy, provided that their length is much 
greater than the cross section diameter but smaller than 
exchange length, this phenomenon is governed by ther- 
mal nucleation and propagation of soliton-antisoliton 
pairs; the associated characteristic time is expected to 
follow an Arrhenius law**'^. For genuine one dimensional 
(ID) Ising systems with single spin- flip stochastic dy- 
namics, a slow relaxation of the magnetization was first 
predicted by Glauber— in 1963. Through Glauber's ap- 
proach, many physical systems were investigated, ranging 
from dielectricgii^i^ to polymersiiiSiii. More fundamen- 
tally, this model has been employed to justify the use 
of the Kohlrauch- Williams- Watts functio n^°i^^ (stretched 
exponential) to fit the relaxation of generalized ID spin 
systems. Also the universality issue of the dynamic crit- 
ical exponenfl^iiiii^ii^iii of the Ising model^^, as well as 
strongly out-of-cquilibrium processes (magnetization re- 
versals^, facilitated dynamicsiS, etc.) have been studied 
moving from the basic ideas proposed by Glauber. 

In this paper, single spin-flip Glauber dynamics is 
used to investigate theoretically the slow dynamics of 
the magnetization reversal in molecular magnetic sys- 
tems. In particular, we extend Glauber's theory^ to 
the Ising coUinear ferrimagnetic chain, as well as to the 
case of a chain in which reciprocal non-collinearity of lo- 
cal anisotropy axes and the crystallographic (laboratory) 
frame is encountered. Such extensions are motivated by 
the fact that (i) in molecular-based realizations of SCM's, 
antiferromagnetic coupling typically has a larger inten- 
sity than the ferromagnetic one; in fact, the overlapping 



of magnetic orbitals, which implies antiferromagnetic ex- 
change interaction between neighboring spins, can be 
more easily obtained than the orthogonality condition, 
leading to ferromagnetisir^2i2ii^; (ii) non-collinearity be- 
tween local anisotropy axes and the crystallographic (lab- 
oratory) frame takes place quite often in molecular spin 
chains. Besides magnetization reversal, the dynamic re- 
sponse of the system to a weak magnetic field, oscillating 
in time at frequency uj, is also studied. Depending on the 
specific experimental geometry, selection rules are found 
for the occurrence of resonant behavior of the a.c. sus- 
ceptibility as a function of temperature (stochastic reso- 
nance) at low frequencies, as well as for slow relaxation 
of the magnetization in zero field at low temperatures. 

The paper is organized as follows. In Sect. II we ex- 
tend Glauber's theory^, originally formulated for a chain 
of identical and coUinear spins, to the more general model 
of a chain with non-collinear spins, possibly with Lande 
factors that vary from site to site. In Sect. Ill we use 
two different theoretical methods (the Generating Func- 
tions approach and the Fourier Transform approach) to 
investigate the relaxation of the magnetization after re- 
moval of an external static magnetic field, starting from 
two different initial conditions: fully saturated or par- 
tially saturated. In Sect. IV we calculate, in a linear ap- 
proximation, the magnetic response of the system to an 
oscillating magnetic field. For a chain of N spins, the a.c. 
susceptibility is expressed as the superposition of N con- 
tributions, each characterized by its time scale; through a 
few simple examples, we show that, depending on the ge- 
ometry of the system {i.e., the relative orientations of the 
local easy anisotropy axes and of the applied field), dif- 
ferent time scales can be selected, possibly giving rise, for 
low frequencies, to a resonant peak in the temperature- 
dependence of the complex magnetic susceptibility. In 
Sect. V we show that the theory provides a satisfactory 
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account for the SCM behavior experimentally observed 
in some magnetic molecular chain compounds, charac- 
terized by dominant antiferromagnetic exchange interac- 
tions and non-coUinearity between spins. Finally, in Sect. 
VI, the conclusions are drawn and possible forthcoming 
applications are also discussed. 



II. THE NON-COLLINEAR ISING-GLAUBER 
MODEL 

In a celebrated paper—, Glauber introduced, in the 
usual ID Ising modeli^, a stochastic dependence on the 
time variable t: i.e., the state of a spin lying on the k- 
th lattice site was represented by a two- valued stochastic 
function (Tfe(i) 

N 

= - ^ ( J/afedfe+i 5MsHe~'"Vfc) , ak{t) = ±1. 
fc=i 

(1) 

J I is the exchange coupling constant, that favors nearest 
neighboring spins to lie parallel (J/ > 0, ferromagnetic 
exchange) or antiparallel (J/ < 0, antiferromagnetic ex- 
change); g is the Lande factor of each spin, and /is the 
Bohr magneton. In the original paper— a ID lattice of 
equivalent and coUinear spins was studied; there the re- 
sponse to a time-dependent magnetic field H(t), applied 
parallel to the axis of spin quantization and oscillating 
with frequency lu, as in typical a.c. susceptibility exper- 
iments, was also considered. 

In order to investigate the phenomena of slow relax- 
ation (for H=0) and resonant behavior of the a.c. sus- 
ceptibility (for 0) in molecular SCM's, we are going 
to generalize the Glauber model (H]) accounting for non- 
collinearity of local anisotropy axes and crystallographic 
(laboratory) frame. To this aim, we adopt the following 
model Hamiltonian 



N 



(Tk(Tk+i + Gk^J.B'^le CTfc) , (Tfe(t) = ±1 



k=l 



(2) 

Jj is an effective Ising exchange coupling that can ap- 
proximately be related to the Hamiltonian parameters 
of a real SC M^'^i^^ : see later on the discussion in Sec- 
tion V. Like in the usual Ising-Glauber coUinear model 
dl]), the spins in Eq. ([2]) are described by classical, one- 
component vectors that are allowed to take two integer 
values (Jfe(i) = ±1, but now the magnetic moments may 
be oriented along different directions, ik, varying from 
site to site. Within this scheme, the Lande tensor of a 
spin on the k-th lattice site has just a non-zero compo- 
nent, gl, along the local easy anisotropy direction z^. 
Denoting by en the direction of the oscilla.ting magnetic 
field, H(i:)=IIe~*'^*eH, we define the generalized Lande 
factor Gfc appearing in Eq. ^ as 



Gk 



9k z/c • en 



(3) 



i.e., a scalar quantity that varies from site to site. Fol- 
lowing Glauber^, we define the single spin expectation 
value Sk{t) — (o'k)t, where brackets denote a proper en- 
semble average, and the stochastic magnetization along 
the direction of the applied field 



N 



N 



^GfcSfc(t). (4) 



k=l 



k=l 



The basic equation of motion of the Glauber mode l^^i^^ 
reads 



-2 (o-fcWcr.^-o 



(5) 



in which w 



represents the probability per unit 
time to reverse the /c-th spin, through the flip -t-Ufc — > 
~(Jk. For a system of N coupled spins, this probability 
is affected by the interaction with the other spins, with 
the thermal bath and, possibly, with an external mag- 
netic field. Among all possible assumptions for the tran- 
sition probability w^i^^-ak as a function of the + 1 
variablesiiii^i^^i^^ { ai , . . . , dfe , . . . , ctat , i } , again following 
Glauber— we require Wak^-ak to be independent of time 
and to depend only on the configuration of the two near- 
est neighbors of the fc-th spin. In zero field, these re- 
quirements are fulfilled by 



w, 



.,H=0 



1 



1 - -jio'k {ak-1 + (Tk+i) 



while in the presence of an external field 



w 



H=0 



(1 - 5k(7k) 



(6) 



(7) 



is usually chosen; the attempt frequency (i.e., the 
probability per unit time to reverse an isolated spin) re- 
mains an undetermined parameter of the model; 7 ac- 
counts for the effect of the nearest neighbors; the param- 
eters 5k have the role of stabilizing the configuration in 
which the k-ih spin is parallel to the field, and destabiliz- 
ing the antiparallel configuration. Thanks to the partic- 
ular choices (0) and ([7]) for the transition probability, by 
imposing the Detailed Balance conditions^ it is possible 
to express 7 and 5k as functions of the parameters in the 
spin Hamiltonian ^ 



7 = tanh {2[3Ji) , 4 ^ tanh {(iGk^lB^) 



(8) 



where /3 = -j^Kp is the inverse temperature in units of 
Boltzmann's constant. Another advantage of Glauber's 
choices ^ and ^ is that the equation of motion ^ 
takes a simple form. In particular, for H = 0, Eq. ^ 
with the choice ([6]) becomes 



ds{t) 
dt 



-aAs{t) , 



(9) 



where s(t) denotes the vector of single-spin expectation 
values {si(t), S2{t),--- , SAr(t)} and A is a square A^ x A^ 
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symmetric matrix, whose non-zero elements are Afc^fc=l 
and Afc_fc_i=Afe_fe+i=--^, with Ai_jv=Ajv,i=-f if peri- 
odic boundary conditions are assumed for the iV-spin 
chain. A closed solution of this set of first-order differen- 
tial equations can be obtained expressing the expectation 
value of each spin, Sk(t), in terms of its spatial Fourier 
Transform (FT) 



Skit) 



— iqk —Xnt 
SqC ^ e ^ . 



(10) 



Substituting Eq. (fTO]) into Eq. ([9]), one readily obtains 
the dispersion relation 



Xq ^ a{l ~ jcosq) , q 



27r 



(11) 



with n — 0, 1,..., N — 1^^. For ferromagnetic cou- 
pling (J/ > 0, hence 7 > 0) the smallest eigenvalue 
Ag=o — ct(l ~ 7) occurs for n=0, independently of the 
number of spins in the chain. For antiferromagnetic 
coupling (J/ < 0, hence 7 < 0) and N even, the smallest 
eigenvalue Ag^jr — ck(l — I7I) occurs for n=^; while in 
the case of N odd, the smallest eigenvalue corresponds 
to a [1 — I7I cos (■^)] , thus depending on the number of 
spins in the antiferromagnetic chain^'*. The characteris- 
tic time scales of the system, r^, are given by 



1 



a{l — 7 cos q) 



(12) 



At finite temperatures the characteristic times Tq are fi- 
nite because I7I < 1; for T — > one has that 1 — I7I van- 
ishes irrespectively of the sign of J/, because 7 jjjj — 
±1. Thus, for H= 0, there is one diverging time scale 
in the T — > limit: t^^o for ferromagnetic coupling and 
Tg=^ for antiferromagnetic coupling (and even N). In the 
presence of a non-zero, oscillating field H(i) = He~'"*, 
the equation of motion (O with the choice ([7]) takes a 
form (see Eq. [M] in Section IV later on) which can still 
be solved, though in an approximate way^, for a suffi- 
ciently weak intensity of the applied magnetic field. 



III. RELAXATION OF THE MAGNETIZATION 
IN ZERO FIELD 

The original Glauber model was formulated for a chain 
of coUinear spins with the same Lande factors: i.e., in 
Eq. ^ one has Gk — g, — 1, ■ ■ ■ , N . Assuming 
that the system has been fully magnetized by means of 
a strong external field, one can study how the system 
evolves if the field is removed abruptly. This corresponds 
to take a fully saturated initial condition 



Sfc(O) = 1, Vfc. 



(13) 



In ferrimagnetic chains, on the other hand, a "par- 
tial" saturation can be reached, provided the antiferro- 
magnetic coupling (J/ < 0) between nearest neighbors is 



"strong enough" , in a sense that will be clarified later on. 
In fact, if the Lande factors for odd and even sites are not 
equal (go ge), through the application of an opportune 
field the sample can be prepared in a configuration with 



Sfc(O) = +1, for fc = 2r + 1 (fc odd) 
s/c(0) = —1, for k ~2r (fc even). 



(14) 



With respect to the case considered by Glauber, it is 
convenient to separate explicitly the expectation values 
of the odd sites, S2r+i{t), from those of the even sites, 
S2r{t). Thus, for H = 0, the set of N equations of motion 
^ can be rewritten as 



TtS2r+l = -a [S2r+1 



a [S2r + 57(s2r+l + S2r-l)] 

^l{s2r + S2r-2)] 



(15) 



In the following, the solutions of (fT5|) will be found using 
two different approaches that yield identical results. 



A. The Generating Function approach 

The Generating Function approach, which closely fol- 
lows the original Glauber's paper, is exposed in detail in 
Appendix. Here, in order to distinguish between the fer- 
romagnetic and ferrimagnetic relaxations, we specialize 
the general solution, Eq. (|A5|1 and Eq. (jA6p . to the two 
different kinds of initial conditions, Eq. (|13p and Eq. (I14p . 
In both cases, we will assume that the exchange cou- 
pling Jj is negative. The "partially saturated" configu- 
ration, Eq. (jl4p . reflects a typical experimental situation, 
in which the antiferromagnetic coupling is much bigger 
(J/ « 100-^1000 K) than the Zeeman energy associated 
with accessible magnetic fields. On the other hand, the 
initial configuration with all the spins aligned in the same 
direction, Eq. (|13p. clearly reflects the experimental situ- 
ation of a fully saturated sample. This condition is easily 
obtained for ferromagnetic coupling (J/ > 0), while it 
may require very strong fields (eventually unaccessible) 
for antiferromagnetic coupling (J/ < 0). 

Let us start from the saturated configuration, Eq. p3|l . 
Substituting the initial condition Sk{0) = 1 for all fc in 
both (E5l) and (M\i. we obtain 



S2r(<) = e "* J2 [^2(r-m)iiat) +l2(r~m)~liiat)] 
m— — 00 

+ 00 



Hence, exploiting the property (|A4p of the Bessel func- 
tions (taking y = 1), and redefining the sums by a unique 
index j , we get 



+00 



-a(l-7)t 



'=~Too (16) 

j=-oo 
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This means that, starting with all the spins aligned in the 
same direction, each spin expectation value (both on even 
and odd sites) decays obeying a mono-exponential law 
with relaxation time t^^q = [q^(1 ~ 7)]~ : which is just 
the characteristic time scale obtained as the inverse of 
the dispersion relation Xq with zero wave number q — 0, 
see Eq. (fTTj) . Notice that t^^o can diverge for T ^ only 
in the case of ferromagnetic coupling, J/ > (7 > 0). 

Let us now consider the partially saturated configu- 
ration, Eq. HH), in which Sfe(O) = 1 for fc odd and 
•Sfc(O) = — 1 for k even. Substituting these initial con- 
ditions in both (jASP and (|A6ll . we obtain 

{+00 
m— — 00 
+00 
S2r+l(i) = e""' [^2(r-m)(7a*) -2'2(r-m) + l(7aO] 

m— — 00 

and, still exploiting the property (jA4p (but now for y = 
— 1), we get 

{-I-CX3 
S2r(t) = -e-"* E (-1)%(7«^) = -e""^'+^^' 
j = -co 
+00 
S2r+i(t) = e-"* E (-lpZ,(7ai)-e-"(i+'^)*. 
j=-oo 

(17) 

Also in this case all the spins of the system relax with 
a mono-exponential law, but now the relaxation time is 
Tq^TT = [«(1 + 7)]~"^j which corresponds to the inverse of 
the eigenvalue Xq with wave number q = n, see Eq. 
Notice that r^^Tr can diverge for T ^ only in the case 
of antiferromagnetic coupling, J/ < (7 < 0). 

Summarizing, according to the sign of the exchange 
constant, both time scales Tq^o (for Jj > 0) and t^^tt 
(for Jj < 0) diverge in the low temperature limit T ^ 0, 
following an Arrhenius law 

r = ^e^^l^^l (18) 

with energy barrier 4|J/| (slow relaxing mode). It is 
worth noting that the remaining relaxation times, given 
by the inverse of the eigenvalues in Eq. (fTTj) with q ^ 
and q TT, always remain of the same order of mag- 
nitude as a"^ (fast relaxing modes). This time scale is 
typically very small {'^ ps) in real systems, and negligible 
with respect to the characteristic times involved in any 
experimental measurement we refer to. 

B. The Fourier Transform approach 

The solutions, and ^T7\ . to the set of equations 
(|15p can alternatively be deduced within the Fourier 
Transform (FT) formalism, which has already been ex- 
ploited to obtain the dispersion relation (fTTj) . Recalling 
the definition ^TU\i of Sk{t) and its spatial FT 

k 



we evaluate Sq at time t = 0, Sq = jjJ2k , 
for the two initial conditions of interest, flS]) and (|14p . 
Starting from the all-spin-up configuration, Eq. (fT3|) . we 
have 

k 

Hence the solution for the expectation value of a spin 
localized on the k lattice site at time t is 

Skit) - Koe'^'e-^-' = e-^'=° * , (21) 

which is identical to since Xq=o — a{l — j). 

Starting from the partially saturated configuration, 
Eq. p4)) . it is useful to rewrite it as Sfc(O) — -e"^ , so 
that the FT at t = is 

fc 

Hence the solution is readily obtained 

Sfc(t) = - ^ 55,^e^9fee-^<.* = _e^^fce-^«='*, (23) 
9 

which is identical to p7|) since A^^jr = ol(\ +7). 

Finally, we observe that Eqs. (PTjl and (I^H)) hold even 
for a ring with a finite number N of spins, while Eqs. 
(I16p and p7p were obtained in the infinite-chain limit. 

C. Slow versus fast relaxation of the spontaneous 
magnetization 

The expectation values, Sfc(i), of spins localized on 
the even and odd sites of a linear lattice at time i, 
computed either with the Generating Functions or the 
Fourier Transform approach, have been shown to dis- 
play a mono-exponential relaxation, see Eqs. p6p and 
(fT7)) . with different time scales, t^^q — [«(1 ^ 7)] "'^ a-nd 
T^=7r[Q!(l + 7)]"^ respectively, depending on the differ- 
ent initial conditions, Eqs. and (dH). As a conse- 
quence, also the macroscopic magnetization, expressed 
by Eq. displays the same mono-exponential relax- 
ation as the single site quantities Sfc(i). 

A chain in which all the magnetic moments are equal 
can be prepared only in the saturated initial configu- 
ration, Eq. (fT5|) . with all the spin aligned in the same 
direction, through the application of an external field. 
Thus, when the field is abruptly removed, such a sys- 
tem will relax slowly at low temperature only if the ex- 
change coupling is ferromagnetic (J/ > 0). In contrast, if 
the exchange coupling is antiferromagnetic (J/ < 0) and 
the chain is "forced" in the saturated state by a strong 
applied magnetic field, the system will relax very fast 
(in a typical time of the order of a~^) when the field 
is removed. Let us now discuss how these results, first 
obtained by Glauber—, are generalized to the case of a 
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chain in which the magnetic moments are colhnear, but 
not equal on each site. 

As pointed out in the introduction, a model with an- 
tiferromagnetic coupling (J/ < 0) but non-compensated 
magnetic moments on the two sublattices is more akin 
to real SCM'sii^ii^S. Yet it is very interesting since, de- 
pending on the intensity of the applied magnetic field, 
the system can be prepared either in the saturated ini- 
tial configuration, Eq. p3p . where all spins are parallel 
to each other, or in the partially saturated one Eq. ([H]), 
where nearest neighbors are antiparallel. In the former 
case, a very strong field is required in order to overcome 
the antiferromagnetic coupling between nearest neigh- 
bors; once the field is removed, the relaxation of the 
magnetization is expected to be fast at low temperatures, 
on the basis of the solution pB]) . In the latter case, the 
partially saturated initial configuration can easily be 
obtained through the application of a smaller, experimen- 
tally accessible magnetic field; when the field is abruptly 
removed, the relaxation is expected to be slow accord- 
ing to the solution (fT7)) . The solution ([T7)l justifies the 
observation of SCM behavior in ferriniagnetic quasi-lD 
compounds like CoPhOMe^S (see Sect. V). 

Summarizing, we have found that when a colhnear 
ferrimagnetic chain is prepared in an initial state - 
fully or partially saturated depending on the intensity 
of the applied magnetic field - once the field is re- 
moved abruptly, the spin system can show fast or slow 
relaxation, respectively. Fast relaxation corresponds to 
stronger fields; unfortunately for the quasi-lD chain com- 
pound CoPhOMe22ii2£, the antiferromagnetic exchange 
constant is so large (|J/| ^ lOOK) that the realization 
of a fully saturated initial configuration would require a 
very high, almost unaccessible field (^ 1000 kOe). Thus 
this compound is not a good candidate for such a kind 
of experiments^. 



IV. MAGNETIC RESPONSE TO AN 
OSCILLATING MAGNETIC FIELD 

In the presence of a magnetic field H, the transition 
probability to be put in the equation of motion is 
w^^^_^^, defined in Eq. ([7]). One obtains 

+ — [(o-feO-fc+i)f + (o-fc-icrfc>J - (5fc| (24) 

that differs from Eq. considered earlier for 11= 0, in 
the presence of both a non- homogeneous term, dk, and 
the time-dependent pair-correlation functions (ct/cCT/c±i)j. 
The latter ones, assuming that the field is so weak to 
induce just small departures from equilibrium, can be 
approximated by their time-independent counterparts^ 



As it is usual in a.c. susceptibility measurements, 
we also assume the time-dependent magnetic field 
H(/;)=IIe^"^*eH, oscillating at frequency w, to be weak 
so that the Sk parameters can be linearized 



tanh(/3/iBGfcH(i)) «/3/isGfcH(t). 



(26) 



The system of equations of motion (I24|) then takes the 
form 



dskjt) 
d{at) 



-Sk{t) + - [Sk+l{t) + Sk-l{t)] 

f3f{pJI)^iBGkii{t) , 



where 



1 -r?^ 
1-1-77^ 



(27) 



(28) 



is a function of the reduced coupling constant pjj and 
we have taken into account that 7 = 2ry/(l + 77^). After 
a brief transient period, the system will reach the sta- 
tionary condition in which the magnetic moment of each 
spin oscillates coherently with the forcing term at the 
frequency co. Expressing the expectation value of a spin 
on the k-th lattice site, Sfe(i), through its spatial FT, Sg, 
the trial solution is 



.Sfc(t) = ^V''e-'"*, 
q 



Substituting the latter in the system ([27|) we get 

aG 

= /3/(/3J/)mbH— 

a[l — 7 cos q) — lU! 



where Gg is the FT of Gfe: 



1 ^ 



(29) 



(30) 



(31) 



The average of stochastic magnetization can readily be 
obtained from ^ as 

N 

(Af), = fiBe~'^' E GgVe^^'e*"''' > (32) 

k—l qq' 

which accounts for non-coUincarity of local anisotropy 
axes with respect to the field direction. Performing the 
sum over all the lattice sites [k indices) yields a factor 
N5q-qi in Eq. ([5^ : substituting the expression ([50]) for 
Sq, one obtains 



GqGq'Sq^-qi 



qq 



■f a{l — 7 cos q) — ioj 



(33) 

Then, considering that GqG^q — |Gg|^, the a.c. suscep- 
tibility is finally obtained dividing ([55)1 by He"*"* 



((Tfc(Tfc+i}t = {ak-icrk)t ~ tanh(/3J/) = r/. 



(25) 



/(a;,r)=iV^|/3/(/3Jz)^— ■ 



(1 — 7 cos q) — iui 



(34) 
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FIG. 1: Temperature dependence of the imaginary part of the complex susceptibility, Eq. (|38[) . of a collinear one-dimensional 
Ising model with alternating spins. Resonant behavior in response to an oscillating magnetic field is possible, at low frequency, 
only when magnetic moments are uncompensated (a,c,d), while a broad peak is found when the net magnetization is zero (b). 
(The curves refer to reduced frequency uj/a = 0.001) 



In principle, the a.c. susceptibility of a chain with N 
spins admits N poles, corresponding to the N eigenvalues 
Xq in Eq. (fTTj) . Each mode is related to a different time 
scale Tg = 1/Xq. In practice, not all the time scales will be 
involved in the complex susceptibility x(a;,r), but only 
the ones selected by Gq. A result similar, at first glance, 
to Eq. (l34|) was deduced by Suzuki and Kubo^*', but in 
their case the relationship was between the time scale Tq 
and the wave- vector-dependent susceptibility x((7,w). In 
contrast, in an a.c. susceptibility experiment only the 
zero-wave- vector susceptibility xi<} — Oj ^) is accessible; 
the peculiarity of Eq. is that other time scales, dif- 
ferent from Tq^o, can be selected thanks to the depen- 
dence of the gyromagnetic factors Gk and of the local 
anisotropy axes on the site position k. This is the main 
result of our study and will be clarified hereafter through 
a few examples. 



A. The a.c. susceptibility of a collinear Ising 
ferrimagnetic chain 

Let us start considering the case of a one-dimensional 
Ising model with two kinds of spins (aligned parallel on 
antiparallel to the chain axis) alternating on the odd 
and even magnetic sites of the lattice with Lande fac- 



tors G2r+i = go and G2r = ffe (integer r), respectively. 
Strictly speaking, a collinear Ising ferrimagnet is char- 
acterized by an antiferromagnetic coupling J/ < 0, but 
also the case J/ > can be treated through Eq. (|34p . 
In fact, since the local axis of anisotropy has the same 
direction for all the spins, the FT of the site-dependent 
Lande factor is 



N/2 



(5e 



N/2 



-iq2r 



(35) 



Taking into account that, in the presence of periodic 
boundary conditions, one has 



N/2 



q2r 



N 



(Sqfl + 5q^^), 



(36) 



it follows that the only non-zero values of Gq are for q = 
and q = TT 



^l = \ [(ffe + 50)^9,0 + {ge - go)Sq^„] 



(37) 



Thus, according to Eq. the parallel a.c. suscepti- 

bility (II = zz) of a collinear Ising chain with alternating 



7 



spins IS 



n-2 



X\\i^,T) 



NfilPfiPJi)^ 
(5e + go)' 



(.9e - .9o)^ 



(l-7)-*(|) (l + 7)-*(^) 



(38) 



It appears that both the relaxation times obtained by 
Suzuki and Kubo^i for the ordinary and the staggered 
susceptibihty of the usual Ising model, namely Tq=o = 
[a{l — 7)]^^ and r^^^ = [a(l + 7)]""^ respectively, do 
coexist in the a.c. susceptibility ((38)) . Notice that, in 
the cj — > limit, the static susceptibility of the Ising 
ferrimagnet in zero fieldii is recovered from Eq. ([35]), 
since one has = i-I^^ = e±2/3./r. 

As regards the dynamic response of the system to an 
oscillating magnetic field applied along the chain axis, 
depending on the sign of the effective exchange coupling 
constant J/, the ferromagnetic {ge + go) or the antifer- 
romagnetic {ge — go) branch of the parallel susceptibility 
([38)) are characterized by a diverging time scale at low 
temperature. In particular, for a coUinear Ising ferri- 
magnet one has J/ < 0, so that r^^^r is diverging, while 
T^=o is short (of the order of a^^, the attempt frequency 
of an isolated spin). Thus, for J/ < 0, a resonant be- 
havior of the a.c. susceptibility versus temperature (at 
low frequencies u>/a <C 1) can only be observed in the 
case ge ^ go (see Fig. Id) when magnetic moments are 
uncompensated, while a broad peak is found in the case 
ge = go when the net magnetization is zero (see Fig. lb). 
Clearly, for J/ > 0, a resonant peak is found in both 
cases (see Fig. la and Ic), because a net magnetization 
is always present in the system. 

Such a resonant behavior of the a. c. susceptibility ver- 
sus T, in ferromagnetic'^^ as well as in ferrimagnetioii 
Ising chains with single spin-flip Glauber dynamics, is a 
manifestation of the stochastic resonance phenomenon^: 
i.e., the response of a set of coupled bistable systems to a 
periodic drive is enhanced in the presence of a stochastic 
noise when a matching occurs between the fluctuation- 
induced switching rate of the system and the forcing fre- 
quency. In a magnetic chain, the role of stochastic noise 
is played by thermal fluctuations and the resonant peak 
in the temperature-dependence of the a.c. susceptibil- 
ity occurs when the statistical time scale, associated to 
the slow decay of the magnetization, matches with the 
deterministic time scale of the applied magnetic field 



^qi.'-^peak) ~ 

UJ 



(39) 




n=3 



n=6 



> 




FIG. 2: Thick arrows denote the projections on the xy plane, 
perpendicular to the chain (helix) axis z, of magnetic mo- 
ments in a one-dimensional Ising helimagnet, for different fold 
symmetries (n — 2,3,4,6). Dashed lines are the projections 
of the local axes of anisotropy, . 



spins experience different fields because of the geometri- 
cal arrangement of magnetic moments. In the following 
we will make the approximation that the Lande tensor 
of a spin on the fc-th lattice site has just a non-zero com- 
ponent g along the easy anisotropy direction z^, so that 
Gk — gzk-eu (see Eq. ([3])). In the crystallographic frame 
{x,y, z), the directors Zk read (integer k) 



( 



27rfc 



\ n 



h cos 6ez. 
(40) 

Let us first consider the case of an oscillating magnetic 
field H applied parallel to z, the helix axis. All the spins 
actually undergo the same field, and since Gk = g cos 9 
independently of the lattice site k, the only peak in the 
FT Gq occurs at q — 



Gl = 5 cos 6* (5g,o (V n). 



(41) 



Following the same procedure as in the previous para- 
graph, the parallel a.c. susceptibility (|| = zz) takes an 
expression (valid for any value of the fold index n of the 
helix) 



cos^ ' 



X„(c.,r) = N^x%f3 /(/3J,)^-^-^_^ (Vn) 



(42) 



B. The a.c. susceptibility of an n-fold helix 

Next, as an example of a non-coUinear spin arrange- 
ment, we consider a system of spins with the local axes 
of anisotropy arranged on an n-fold helix (see Fig. 2); 6 
is the angle that the local axes form with z, the unique 
axis of the helix (i.e., the chain axis). In this case the 
Lande factors are equal on all lattice sites, but different 



that differs from Glauber's result for the coUinear Ising 
chain- only by the geometrical factor cos^ 9. For fer- 
romagnetic coupling, Ji > 0, the relaxation time tq = 
[a(l — 7)]"^ diverges as T ^ 0, and a resonant behavior 
of the a.c. parallel susceptibility versus temperature is 
found, at low frequency, when the oscillating field is ap- 
plied parallel to the helix axis, z, along which spins are 
uncompensated: see Fig. 3a, which refers to the case of a 
two-fold helix [n = 2) . 
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FIG. 3: (color online) Temperature dependence of the imaginary part of the parallel (|42p and perpendicular (|45p complex 
susceptibility of an Ising chain with two-fold helical spin arrangement. The local axes zi and Z2 were assumed to form an angle 
= with z, the chain axis (unique axis of the helix). Different curves refer to different values of u/a: 0.0001 (continuous, red 
line); 0.0002 (dashed, green line); 0.0005 (dashed single-dotted, blue line); 0.0010 (dashed double-dotted, violet line). Resonant 
behavior in response to an oscillating magnetic field is possible, at low frequency, only for field applied in a direction where 
magnetic moments are uncompensated (a,c), while a broad peak is found (b,d) when there is no net magnetization along the 
field direction. 



Let us now consider the case of an oscillating magnetic 
field H applied perpendicularly to the chain axis. In this 
configuration, it is useful to distinguish the case n = 2 
from the general case n > 2. 

• 71^2 

In this case, it is vi^orth noticing that for H parallel 
to y, one has identically Gr = for any lattice 
site r. Thus, = and the corresponding a.c. 
susceptibility is identically zero 



Xyy{uJ,T) = 



(43) 



(not shown). In contrast, for H parallel to x, one 
has Gr = ^gsinO on odd sites and Gr — +gaiTi9 
on even sites. The FT is 



^ N/2 



_g-ig(2r-l) _|_ g-j«2r\ 



gsin6'-((5,^o + <5g,^)(l 



(44) 



where we have taken into account Eq. ([55]) . Thus, 
for ferromagnetic coupling, J/ > 0, the relaxation 



time — [a(l + 7)] ^ does not diverge as T 
and the perpendicular a.c. susceptibility 



0, 



X.A^,T)^Nfil(3f{pji) 



g"^ sin^ I 



(l+7)-»(S) 



(45) 



does not present a resonant behavior as a function 
of temperature; rather, it presents a broad max- 
imum (see Fig. 3b). Clearly, in the case of an- 
tiferromagnetic coupling, J/ < 0, the behavior of 
the susceptibility components is reversed: a broad 
maximum is found for the temperature dependence 
of the parallel susceptibihty Xzz (w, T) (see Fig. 3d), 
while a resonant behavior is found for the perpen- 
dicular susceptibility Xxx{^,T) (see Fig. 3c). 

• n> 2 

In this case, denoting by e^; • en and e.y ■ en the 
directors of the in-plane field, the FT's of Gk are 
given by 



Gy 



= (ex • eH)gsin6'i {5^ 2^ 
^^(e,.eH)gsin0E£iSin(^)e 
= (gj, • eH)gsin6'Jj {5^ 21 



-iqk 



-iqk 



(46) 
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Remarkably, as just \Gq\'^ appears in Eq. p4|) . the 
general result for the in-plane susceptibility turns 
out to be independent of the field direction. Thus, 
for n> 2, the perpendicular (_L) a.c. susceptibility 
of the n-fold helix is given by 

1 

X -SlYl^e ^75-^ , (47) 

2 [l-^cos(2^)]_z(H)' ^ ^ 

where we have exploited the fact that cos ^) = 
cos (^) for the term appearing at the denominator 
of Eq" 

Summarizing, in the general case of an Ising chain 
with an n-fold helical spin arrangement {n > 2), we 
have explicitly shown that a resonant behavior of the 
a.c. susceptibility versus temperature, similar to the one 
displayed by ferromagnetio^*^ and ferrimagnetioii Ising 
chains with coUinear spins, is possible only for field ap- 
plied in a direction where magnetic moments are uncom- 
pensated. In contrast, a broad peak is found when there 
is no net magnetization along the field direction. 

V. APPLICATION TO REAL SINGLE CHAIN 
MAGNETS 

In this Section we will apply the developed formalism 
to some real compounds as representative realizations of 
SCM's; for the three selected systems - we know this 
restriction is far from being exhaustivei*^ - a.c. suscep- 
tibility data on single crystal are available, which is a 
fundamental requirement for checking the proposed selec- 
tion rules. The considered system o^^'^^i'^^ are character- 
ized by the alternation of two types of magnetic centers 
along the chain axis, so that at least two spins per cell 
have to be considered; moreover, the magnetic moments 
are not coUinear, the dominant exchange interactions are 
antiferromagnctic and a strong single-ion anisotropy is 
present, which favors magnetization alignment along cer- 
tain crystallographic directions Zk ■ The static properties 
of these compounds, like magnetization and static sus- 
ceptibility, are generally well described using a classical 
Heisenberg model with an isotropic exchange coupling J 
and a single-ion anisotropy D. Thus, in order to describe 
the dynamic behavior in response to a weak, oscillating 
magnetic field by means of the previously developed the- 
ory, it is necessary to relate the Hamiltonian parame- 
ters of such a classical spin model to the exchange con- 
stant J] of the effective Ising model ([2]). In the following 
we will show, through a few examples on real systems, 
that indeed, depending on the geometry, selection rules 
are obeyed for the occurrence of slow relaxation of the 
magnetization at low temperatures (/3|J/| 3> 1), as well 
as for resonant behavior of the a.c. susceptibility as a 
function of temperature at low frequencies. As regards 
the frequencies involved in an a.c. susceptibility exper- 
iment on real SCM's, generallylj"^ they lie in the range 



10""'^ 10^ Hz, while the attempt frequency a is of the 
order of 10^° 10^"^ Hz. Thus, for a typical experiment, 
a resonant peak in the a.c. susceptibility can safely be 
observed provided that at least one of the characteristic 
time scales Tq involved in (j34[) diverges at low T, in order 
for the condition ([55)1 to be satisfied. 

A. The Mn"'-based Single Chain Magnet 

In the one-dimensional molecular magnetic compound 
of formula [Mn(TPP)02PPhH]-H20, obtained by re- 
acting Mn(III) acetate mesotetraphenylporphyrin with 
phenylphosphinic acid'^^-, hereafter denoted by Mn^^^- 
based SCM, the phenylphosphinate anion transmits a 
sizeable antiferromagnctic exchange interaction that, 
combined with the easy axis magnetic anisotropy of the 
Mn^^^ sites, gives rise to a canted antiferromagnctic ar- 
rangement of the spins. The static single-crystal mag- 
netic properties were analyzed in the framework of a clas- 
sical spins Hamiltonian 

N/2 

+ e-^-'^iB^'^g'^P[Si_, + Si]} (48) 

where J < is the antiferromagnctic nearest neighbor 
exchange interaction between S = 2 spins. Z? > is 
the uniaxial anisotropy favoring two different local axes, 
alternating along odd and even sites respectively; both 
axes form an angle 9 = 21.01° with the crystallographic 
c axis, while they form opposite angles of modulus cj) — 
56.55° with the a axis (see Fig. 4). Thus we can write 

Z2r-1 = S\n9 COS(t)Bx ~ SVnO SVUcjiBy + COsOGz and Z2r = 

sin 9 cos 4>ex + sin 9 sin (pey -\- cos 9ez ■ 

A best fit of the static single-crystal magnetic suscep- 
tibilities, calculated via a Monte Carlo simulation^ pro- 
vides J = —1.34 K and D = 4.7 K; the gyromagnetic 
tensor G°'^ is diagonal and isotropic with = 1.97. 
Equivalent results can be obtained calculating the static 
properties of model l[15|) via a transfer matrix approach^!. 
Since the uniaxial anisotropy D is rather strong with re- 
spect to the exchange coupling | J|, as a first approxima- 
tion one can assume the two sublattice magnetizations 
to be directed just along the two easy axes, Z2r-i and 
z;2r, so that the chain system (|48)) can be described by a 
non-coUinear Ising model formally identical to Eq. ([2]), 
with an effective"^^ Ising exchange coupling J/ and a gen- 
eralized Lande factor Gk defined as, respectively 

Jl = JS'(S' + 1) COs(z2r-l • Z2r) 

Gr = gl ^S{S + 1) {ir ■ en). (49) 

Depending on the orientation of the oscillating mag- 
netic field with respect to the crystallographic axes, the 
FT of the generalized Lande factor takes the following 
forms 

Gq = g\\^S{S+l) 
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FIG. 4: (color online) Disposition of local axes (22, -1 and 
Z2r) and magnetic moments (red arrows) in the Mn^'^-based 
real SCM, discussed in Sect. V.A, with antiferromagnetic 
effective Ising exchange coupling J/ < . Right: Schematic 
view of the chain structure {z is the chain axis) along the 
crystallographic s axis. Left: projections of local axes (dashed 
lines) and of magnetic moments (red arrows) in the xy plane, 
perpendicular to the chain axis. 



N/2 



X — ye 
N ^ 



sin 01 cos (/>! 5^,0, H||x 

sin^i sin0i(5g_^, H||y (50) 

cos 01 5^,0, H||z (chain axis) 



The corresponding a.c. susceptibiUty takes the expres- 
sion 



x(u;,T) = 7VM^/3/(/3J,)(gll)2[5(5 



X < 



sin C/sm 



1)] 

(l-7i)-i(^)' II 



(51) 



cos 



(l-7f)-»(^) 



, H||z (chain axis) 



Taking into account that, for the Mn^^^ SCM under study, 
the "true" exchange coupling, J in Eq. (j48|) . is antifer- 
romagnetic, and that the angle between the two easy 
anisotropy axes i\ and Z2 is (5 = 34.6° < 90° (see Fig. 
4, right), from Eq. (|49| it follows that also the effective 
Ising exchange coupling is antiferromagnetic, J/ < 0. As 
a consequence, in the low temperature limit /3| J/j — > cxd, 
the relaxation time r^^^ diverges, while t^=o does not. 
Thus, for low frequencies loja <IC 1, the a.c. susceptibil- 
ity presents a resonant behavior only when the oscillat- 
ing magnetic field is applied along the crystallographic y 



axis, i.e. the direction, perpendicular to the chain axis, 
along which the magnetizations of the two sublattices 
are uncompensated (see Fig. 4). In contrast, when H 
is applied parallel to z (the chain axis) or to .t, namely 
two directions along which the magnetizations of the two 
sublattices are exactly compensated, no resonant behav- 
ior is expected. These theoretical predictions turn out 
to be in excellent agreement with experimental a.c. sus- 
ceptibility data^^ obtained in a single crystal sample of 
[Mn(TPP)02PPhH]-H20, thus confirming that such a 
Mn'^^-based canted antiferromagnet is a bona fide SCM. 



B. The Dy"^-based Single Chain Magnet 

The molecular magnetic compound of formula 
[Dy(hfac)3(NITPhOPh)], hereafter denoted by Dy"i- 
based SCM, belongs to a family of quasi one-dimensional 
magnets in which rare earth ions (with spin S) and or- 
ganic radical ions (with spin s = 1/2) alternate them- 
selves along the chain axis, z, which in this compound co- 
incides with the crystallographic b axis. Static measure- 
ments in single crystal samples suggest^ that there is an 
antiferromagnetic exchange interaction between neigh- 
boring Dy^^^ ions, whose easy anisotropy axes are canted 
with respect to the chain axis in such a way to generate an 
uncompensated moment along b, while the components 
in the ac plane are compensated. Thus, as far as the 
dominant exchange interaction J < between Dy^^^ ions 
is taken into account, the spin Hamiltonian of the sys- 
tem is quite similar to Eq. (|48p . However, with respect to 
the Mn^^^-based chain, the crystal structure of the Dy^^^- 
based SCM is more complicated, not only owing to the 
presence of two kinds of magnetic centers (the Dy'^^^ ions 
and the organic radical ions), but mainly because the 
system is formed by two different families of chains, with 
two almost orthogonal projections of the easy axes in 
the ac plane, perpendicular to the chain axis: this "ac- 
cidental" (in the sense that it is not imposed by symme- 
try) orthogonality is the reason for the nearly isotropic 
magnetic behavior displayed by the system within such 
a plane'^*'. 

We adopt a simplified model formally equivalent to Eq. 
(|48p . Taking into account only the dominant antiferro- 
magnetic exchange interaction (J < 0) between neigh- 
boring Dy^^^ ions (which indeed are next nearest neigh- 
bors in the real system) and their uniaxial anisotropy 
(L* > 0), the system can approximately be described by 
the classical spins Hamiltonian where now l^fcl = 1. 
By means of a classical Transfer Matrix calculation, the 
static properties of the Dy^^^-based SCM turn out to be 
satisfactorily fitted^^ by J = -6 K, D = 40 K, gll = 10, 
with the two easy anisotropy axes Z2r-i, Z2r forming 
equal angles 9 w 75° with the chain axis z. (Notice that 
the latter property holds true for both families of chains.) 
Also in the case of the Dy^^^-based SCM, the uniaxial 
anisotropy D turns out to be sufficiently strong with re- 
spect to the exchange coupling |J| in order to assume, 
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FIG. 5: (color online) Disposition of odd and even local axes 
(z2r-i and Z2r) and magnetic moments (thick arrows) in the 
Dy'^'-based real SCM, discussed in Sect. V.B, with ferro- 
magnetic effective Ising exchange coupling J/ > . Top; 
Schematic view of the chain structure (z is the chain axis), 
displaying the two families of chains (A, with red magnetic 
moments, and B, with green magnetic moments). Bottom: 
projections of magnetic moments in the xy plane, perpendic- 
ular to the chain axis. 



limit /3J/ — > oo, the relaxation time r^^o diverges, while 
Tq=T^ does not. Thus, the a.c. susceptibility 

X(c^,r) « iVM|/3/(/3J,)(gll)2 
„2 a 1 



cos 



sm 



(chain axis) ^^^^ 



(i+7i)-^(^) 



H _L 2 



is expected to have a resonant behavior, for low frequen- 
cies Lo I a 1 , only when the oscillating magnetic field is 
applied parallel to the chain axis, z, along which the mag- 
netizations of the two sublattices are uncompensated (see 
Fig. 5, top). Such a theoretical prediction turns out to 
be in excellent agreement with the experimental a. c. sus- 
ceptibility data^^^ obtained in a single crystal sample of 
[Dy(hfac)3(NITPhOPh)]oo, thus confirming that also the 
Dy^^^-based canted antiferromagnet is a bona fide SCM. 
The only qualitative difference, with respect to the Mn^^^- 
based chain is that, due to the different geometry of the 
spin arrangement and of the local anisotropy axes with 
respect to the chain axis, the resonant behavior of the a. c. 
susceptibility is now observed for field applied parallel to 
the chain axis, rather than perpendicular to it. 



The CoPhOMe (Co"-based) Single Chain 
Magnet 



as a first approximation'^*, the two sublattice magneti- 
zations of Dy^^^ to be directed just along the two easy 
axes. Thus one can define an equivalent non-coUinear 
Ising model where the effective Ising exchange cou- 
pling Jj and the generalized Lande factor Gr are now 
defined as 

J/ = JC0s(z2r-l • Z2r) 

Gr = gl (zr-en). (52) 

Depending on the orientation of the oscillating magnetic 
field with respect to the crystallographic axes, the FT of 
the generalized Lande factor takes the form 

1 

= ^E^~''"^[^''(^2.-l-eH) + (z2.r-eH)] 

II \cos6S„o, Hllz (chain axis) 
a 5" < . . . „ I • (53) 

smf dq^TTi H _L z 

It is important to notice that this result holds true for 
both families (A,B) of chains. Next, we observe that since 
in the Dy^^^-based SCM, the spins on opposite sublattices 
are coplanar with the chain axis, the angle between Z2r-i 
and Z2r is just 26 150° > 90°. Taking into account that 
the "true" exchange constant in Eq. ([48]) is antiferromag- 
netic, J < 0, from Eq. it follows that the effective 
Ising exchange coupling is now ferromagnetic, Jj > (see 
Fig. 5, top). As a consequence, in the low temperature 



In the molecular magnetic compound of for- 
mula [Co(hfac)2NITPhOMe], hereafter denoted by 
CoPhOM o^^'^" , the magnetic contribution is given by 
Cobalt ions, with an Ising character and effective S = 
1/2, and by NITPhOMe organic radical ions, magnet- 
ically isotropic and with s = 1/2. The spins are ar- 
ranged on a helical structure, schematically depicted in 
Fig. 6, right, whose projections in a plane perpendic- 
ular to the helix axis z (coincident with the crystallo- 
graphic c axis), are represented in Fig. 6, left. The prim- 
itive magnetic cell is made up of three Cobalts (black 
arrows) and three organic radicals (red arrows). Al- 
though the effective spins of the two types of magnetic 
centers have the same value, the gyromagnetic factors 
are different: gco 7^ gn', thus, since the nearest neighbor 
(Cobalt-radical) exchange interaction is negative (and 
strong, I J| « 100 K)^, the sublattice magnetizations are 
not compensated along z, whereas they are compensated 
within the xy plane perpendicular to the chain axis z. 
For this compound, which was the first to display SCM 
behavior^^'^*^, static measurements on single-crystal sam- 
ples has not been interpreted in terms of a simple model 
yet, due to the complexity of the system itself. Thus, 
a relationship such as and ([5^ . which associate the 
Ising Hamiltonian ^ parameters (J/ and Gk) with those 
of a more realistic Hamiltonian, is still missing. However, 
the dynamic behavior has been thoroughly investigated 
treating - for the sake of simplicity - both the Co^^ and 
the organic radical spins as Ising variables, with tr = ±1. 
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FIG. 6: (color online) Disposition of even and odd local axes 
(dashed lines) and magnetic moments (thick arrows) in the 
CoPhOMe real SCM, discussed in Sect. V.C, with antiferro- 
magnetic effective Ising exchange coupling J/ < . Right: 
Schematic view of the chain structure {z is the chain axis) 
along the crystallographic y axis. Left: projections of local 
axes (dashed lines) and of magnetic moments (thick arrows) 
in the xy plane, perpendicular to the chain axis. 



The effective Ising Hamiltonian reads 

N/Q 3 
/=1 m=l 

(z2m-eH)]} (55) 

with / magnetic cell index and m site label with boundary 
conditions cr/j — (7/+i^i. Since all the local axes [k — 
1, • • • ,6) form the same angle 9 k, 55° with the z axis, 
when a magnetic field is applied along z, the FT of the 
generalized Lande factor is simply given by 



.9i?.e 



~iq2r 



Gy=cos^l5](gcoe-^(2'-i) 

r— 1 

COS 

= —^[{gco + gR)5qfl + {gco - 9R)^q,^] (56) 



which, except for the prefactor cos 9, is quite similar to 
Eq. p7p for the coUinear Ising chain with alternating 
spins. Thus, the parallel a.c. susceptibility is 



XII (w,T) 



r (3Co+5fl)^ 



cos" 



^(l-7/)-»(^) (l + 7/)-*(^) 



(57) 



Taking into account that the effective exchange coupling 
of CoPhOMe is negative (J/ < 0), the antiferromagnetic 
branch of the parallel susceptibility is characterized by 
a diverging time scale t^^^ = [a(l + 7/)]^^ at low tem- 
perature, so that, for low frequencies uj/a <C 1, X||(<^i2^) 
displays a resonant behavior. 



Let us now consider the case of a field applied in the 
plane perpendicular to the chain axis. For H||x (see Fig. 
6, left) one has, letting = f 



G 



2r-l 



gi^sin6'cos [fco(2r — 1)] 



— 9Co sin 9 cos [fco2r] 
so that the FT takes the form 

1 "^/^ / 
'^9 = sm9-j2e-''^\9cocos{ko2r) 

r=l 

+ [cos(A;o) cos(A;o2r) + sin(A;o) sin(fco2r)] 

= ^sin9[{gco+gRe'^''-''°^)iSg,ko+K^+ko) 
+ (3co + 5fle'(«+'=°))((5,,_fc„ +<5,,^_feJ] 



(58) 



(59) 



where, as usual, we have exploited Eq. (j36p . Thus it 
follows that 



G, 



G. 



4 

sin 



{gco + 9r) 
- (gco - gR)- 



(60) 



The corresponding relaxation times are Tq=±i_ 
and T, 



SO that, summing the four contribu- 



tions we obtain the perpendicular a.c. susceptibility 



{gco + gR)^ , {gco 



sm 



gR? 



(61) 



In conclusion, for the six-fold helix model with alter- 
nating spins and Ising exchange coupling in Eq. (j55p . 
the parallel and perpendicular components of the a.c. 
susceptibility, x\\{^^T) and x^i^^T), display a behavior 
similar to that of a ferrimagnetic chain with alternating 
spins (see (1551) ) and of an n-fold helical spin arrangement 
with equivalent spins (see (|T7|) ). respectively. In spite of 
the approximations involved in model (j55p to describe 
the real CoPhOMe molecular magnetic chain, the two 
calculated susceptibilities (|57|) and ([6T|) , qualitatively re- 
produce the dynamic behavior of this compoun d^^'^*^ . In 
fact, no out-of-phase a.c. susceptibility (imaginary part) 
is observed when the field is applied in the plane per- 
pendicular to the chain axis, z, for the experimental fre- 
quencies (1 10^ Hz)^. In contrast, when the oscillating 
field is parallel to z, a resonant behavior is observed as 
a function of temperature. Even though our theoretical 
treatment holds only for small deviations from equilib- 
rium, it is worth mentioning that the absence of slow 
relaxation for fields applied in the perpendicular plane is 
evidenced in the low temperature magnetization curve as 
well: at low enough temperatures, a finite-area hysteresis 
loop is present only when a static field is applied parallel 
to the chain axis, while no hysteresis is observed in the 
in-plane magnetization curve^Si^. 
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VI. CONCLUSIONS 

In conclusion, in the framework of a one-dimensional 
Ising model with single spin-flip Glauber dynamics, 
taking into account reciprocal non-coUinearity of local 
anisotropy axes and the crystallographic (laboratory) 
frame, we have investigated: (i) the dynamics of mag- 
netization reversal in zero field, and (ii) the response of 
the system to a weak magnetic field, oscillating in time. 
We have shown that SCM behavior is not only a feature 
of coUinear ferro- and ferrimagnetic, but also of canted 
antiferromagnetic chains. In particular, we have found 
that resonant behavior of the a.c. susceptibility versus 
temperature in response to an oscillating magnetic field 
is possible, at low frequency, only for fields applied in a 
direction where magnetic moments are uncompensated. 
In contrast, a broad peak is expected when there is no 
net magnetization along the field direction. 

The role played by geometry in selecting the time 
scales involved in a process is an important and pecu- 
liar result, typical of magneto- molecular approach to low- 
dimensional magnetism. In fact, magnetic centers with 
uniaxial anisotropy usually correspond to building blocks 
with low symmetry22ii^, which - in turn - often crystal- 
lize in more symmetric space groups, realizing a recip- 
rocal non-collinearity between local anisotropy axes as a 
natural consequence^'^. Thus the family of real SCM's, 
to which our model applies, does not restrict to ad-hoc 
synthesized compounds but, instead, is expected to grow 
larger in the future'^. As a validity check of our selec- 
tion rules (as well as a tutorial exemplification), we have 
shown how our theory applies successfully to three differ- 
ent molecular-based spin chains; when possible, we have 
put the parameters of our model Hamiltonian ^ in re- 
lationship with those of more general models, typically 
used to fit the static properties of the corresponding com- 
pounds. Needless to say that the possibility of schema- 
tizing the chosen three compounds with Hamiltonian Q 
relies on the fact that at low enough temperatures they 
behave as chains consisting of two-level units coupled by 
a fully anisotropic exchange interaction. The latter as- 
sumption is expected to hold also for spin larger than 
one-half in the presence of strong single ion anisotropy, 
provided that domain walls still remain sharp^l^. In 
this case each single magnetic center follows a thermally 
activated dynamics, with an energy barrier Aq, and well 
established heuristic arguments''^ suggest to replace the 
attempt frequency by a = aoe~'^'^". 
A naive application of our 3-fold-helix results ([^^ 
and (|47p to the recently synthesized non-coUinear Dya 
cluster would prevent the observation of slow relaxation, 
while Single-Molecule-Magnet dynamics is indeed there 
observed even in the presence of compensated magnetic 
moments^. However, such a behavior in the classical 
regime, i.e. far from level crossings where underbarrier 
processes of quantum origin are important, is observed 
for Dya in non-zero field and the resonant behavior is 
due to a change of the relative population between the 



lowest and the first excited Kramers doublets of each Dy 
ion: For sure this mechanism cannot be accounted for 
when dealing with two-level elementary variables, like ak 
in Hamiltonian Q . An extension of our model to mul- 
tivalued (Tfe's definitely deserve to be considered in the 
next future. 

Beyond molecular spin chains, our approach might also 
be used to model monatomic nanowires showing slow re- 
laxation of the magnetization at low temperatures^ and, 
possibly, one-dimensional spin glasses^ (provided that 
quenched disorder is somehow taken into account). In 
this regard, the question of distinguishing between SCM 
and spin-glass behavior in quasi-lD systems is still a hot 
topic of discussion^i^i^i^. 

After the successful organization of Single-Molecule 
Magnets onto surfaces^^i^Si^, the grafting of properly 
functionalized SCM's on substrates represents a foresee- 
able goal as well as a fundamental step for their possi- 
ble use as magnetic-memory units'^. Technologies em- 
ploying more traditional materials but based on alter- 
native geometrical arrangement of magnetic anisotropy 
axes with respect to the switching field, such as in per- 
pendicular recording^ or processional switching^, are 
already at the stage of forthcoming implementation in 
devices. Were SCM's to be considered as a possible route 
to tackle the main issues of high-density magnetic stor- 
age - i.e. optimization of the signal- noise ratio, thermal 
stability and writability^^ - the proposed selection rules 
for slow relaxation, and related bistability, might find an 
application in magnetic-memory manufacture as well. 
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APPENDIX A: THE GENERATING FUNCTIONS 
APPROACH 

Assuming periodic boundary conditions and defining 
the two generating functions 

-t-oo 

C{y,t) = J2 2/"^+'s2.+i(<) 

r— — oo 
+ 00 

g{y,t) = y'^^'^rit), (Al) 

r— — 00 

Eqs. ([T5)) take the form of two differential equations (with 
the dot indicating the first derivative with respect to the 
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adimensional variable at) 

U{y, t) = -C{y, t) + i7(y + y-^)g{y, t) 
1 Q{y, t) = -Giy, t) + \^[y + y-^)C{y, t) . 



(A2) 



The system (|A2| can be decoupled through the substi- 
tution 



U{y,t) = C{y,t)+g{y,t) 
W{y,t) = C{y,t)~g{y,t), 



(A3) 



from which we directly get 



iuiy,t) = -il-j,piy,t) 
\W{y,t)^-{l + u)W{y,t), 

with ly ~ ^l{y + y^^)- The solutions of these two 
equations are U{y,t) = 0)e-(i-'")"* and Wiy,t) = 
W(j/, 0)e~'^^+'^^"* that, exploiting the property 



+y''-+'giy.0)l2r+i{iat)] 

r— — 00 

+y^^-'£{y,0)l2r-ihat)]. 



According to (|A1[) . now C{y,0) and G{y,0) can be ex- 
pressed again in terms of the initial single spin expecta- 
tion values S2r(0) and S2r+i(0) respectively 



+00 

r— — 00 

+00 

x[2/''' y'"+'s2m+i(0)l2r(7at) 

m—~oo 

4-CXD 



exp 



-f-oo 

E y^Tk{x) (A4) Putting fc' = fc -f to we have 



k— — c 



of the Bessel functions of imaginary argument 1^ (cc) , can 
be rewritten as 



/e— — 00 

W(2/,0 = W(y,0)e-"* ^ (-l)Vlfc(7ai) 

fc^ — c>o 



+00 



+ CX) 



in— — 00 



Performing the inverse transformation of (|A3p . the solu- 
tions to the system (jA2p can be obtained 



1 



Y / [Z^(y,0)Xfe(7at) + (-l)'>V(2/,0)Jfe(7«i)] 

fc — — 00 

+ C)C 

E / [W(y,0)Zfe(7at)-(-l)'=>V(2/,0)Jfe(7«i)] . 

Then, separating the /c-odd from the /c-even terms in both 
sums, we get 

+00 

C{y,t)^e-''' J2 [y^'C{y,Q)l2r(nat) 



[s2m+l{Q)l2(k' ^m){iat) + S2m (0)l2(fc' + l (7"*)] ■ 

Comparing this latter result with the definition of t) 
(jAT|) and requiring for the terms corresponding to the 
same power of y to be equal, an explicit function for the 
odd spin expectation values is readily obtained 

+00 

S2r+l(0 = Y [S2m+l(0)l2(r-m)(7at) 

m— — 00 

+ S2m(0)X2(r-m) + l(7a^)] • (A5) 

Substituting C{y,Q) e g{y,Q) in the solution found for 
Q{y,t) and performing similar passages, we obtain the 
expectation value for even sites 

+00 

S2r(i)=e""* Y [S2m(0)l2(r-m)(7at) 
m— — 00 

+ S2m+l(0)X2(r-m)-l(7ai)] ■ (A6) 
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